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Abstract 

We find the microscopic spectral densities and the spectral correlators associated with multi- 
critical behavior for both hermitian and complex matrix ensembles, and show their universality. 
We conjecture that microscopic spectral densities of Dirac operators in certain theories without 
spontaneous chiral symmetry breaking may belong to these new universality classes. 



CPT-97/P.3526 

NBI-HE-97-45 

NORDITA-97/73 

NSF-ITP-97-137 

hep-th/971200^ 



1 Introduction 



It has been conjectured that the spectra of massless Dirac operators near eigenvalues A = display 
universal features that can be extracted from zero-dimensional theories of random matrices [Q, |2|. The 
natural object to focus on is the so-called microscopic spectral density 

Ps{X) ^ ^^^^p{^) , (1-1) 

which displays the behavior of the spectral density near the origin, as measured on the scale of the 
space-time volume V. 

For QCD-like theories in (2n + l)-dimensions (gauge groups SU{Nc) for Nc > 3, Nf fermions in 
the fundamental representation of the gauge group) the relevant matrix ensemble is that of hermi- 
tian matrices, the unitary ensemble Q. In even space-time dimensions, where chiral symmetry can 
be defined, its possibly spontaneous breakdown is known to be related to a non-vanishing spectral 
density at the origin [^: {^JJ^p) ~ /o(0). In odd space-time dimensions the closest analogue to chiral 
symmetry is parity. This raises the obvious question of whether flavor or flavor-parity symmetries 
can be spontaneously broken in odd dimensions, similar to the breaking of chiral flavor-symmetries 
in even dimensions. While there is yet no direct evidence from, e.g., Monte Carlo simulations that 
this actually occurs for QCD-like theories (see, however, ref. [^), one can entertain this idea, and 
look at some of the consequences. It has been suggested that for an even number of flavors, and 
three space-time dimensions, the pattern of symmetry breaking may be of the flavor-breaking kind 
U{Nf) -^ U{Nf/2) X U{Nf/2) (sj, |2|. An order parameter for such a transition is \T,\ = \{iljip)\, in 
complete analogy with the breaking of chiral symmetries in even dimensions. By the Banks-Casher 
relation 0], this order parameter is proportional to the spectral density at the origin, p{0). In such 
a scenario one can evaluate the microscopic spectral density ps{X) by methods very analogous to the 
even-dimensional case |2|. The result has been found to be & 0] 

PsW = ^vrV(0)2A{j^^.+i(M0)A)2 + J^^._i(^p(0)A)2 

- J^^+i (^p(0)A)J^^_3(7rp(0)A) - J^^,_i (^p(0)A)J^^,+|(7r/,(0)A)} (1.2) 

for the Gaussian distribution. As was recently proven in ref. [0], this microscopic spectral density is 
highly universal. It follows from hermitian random matrix models with arbitrary potentials V{M) that 
support a large- A^ spectral density which is non-vanishing at the origin. It was confirmed in 0] that this 
microscopic spectral density is consistent with the appropriate generalization of the Leutwyler-Smilga 
spectral sum rules [^ to this odd-dimensional situation. 

The question of whether the symmetry breaking U{Nf) -^ U{Nf/2) x U{Nf/2) actually occurs in, say, 
(2+l)-dimensional QCD, is of course a dynamical question that is not addressed in the above consider- 
ations. To achieve such a spontaneous symmetry breaking, the eigenvalues of the Dirac operator must 
accumulate, as the volume V is taken to infinity, sufficiently fast near A = 0. A simple scaling argument 
shows that the accumulation must be such that the average level spacing AA among the eigenvalues 
must be roughly constant near A = 0, and inversely proportional to the volume V. This is precisely 



what motivates the introduction of the microscopic density ps{X) as in eq. (1.1). An obvious question 
to ask, then, is to what extent universal microscopic spectral densities can be extracted if there is 
no spontaneous symmetry breaking. In other words, suppose the dynamics of the massless theory is 
such that the eigenvalues of the Dirac operator do not accumulate sufficiently fast for a condensate 



to form. Then the whole connection to the Leutwyler-Smilga sum rules (and their odd-dimensional 
generalizations) is lost, and the argument in favor of universal microscopic spectral densities has to be 
reconsidered. This question is particularly interesting in the case of (2+l)-dimensions, because there 
the spontaneous breaking of flavor or flavor-parity symmetries is far from obvious. 

But also the even-dimensional case (corresponding to integration over complex matrices, the so-called 
chiral unitary ensemble, |^ for QCD-like theories) merits closer attention in this connection. For 
example, the universal results of refs. [Q, are derived for an arbitrary number of flavors Nj. The 
dependence on Nj is simple, occurring only in the indices of the relevant Bessel functions. This rather 
mild variation as a function of Nf is totally consistent with the fact that the only essential flavor 
dependence enters in the precise chiral symmetry breaking pattern SUL{Nf) ® SUR{Nf) — > SU{Nf), 
and not in detailed dynamical questions such as quark screening etc. Surely the assumption about 
this chiral symmetry breaking pattern must break down as the number of flavors is increased. For 
QCD an almost certain upper limit is (for quarks in the fundamental representation) Nj = 17, for 
which the leading-order coefficient of the beta function has changed sign. The renormalized theory in 
this case is presumable "trivial" , and almost certainly not in a phase of chiral symmetry breaking]^ 
However, although the renormalized coupling for infinite cut-off may be vanishing, it is perhaps not a 
priori obvious that the spectrum of the Dirac operator is identical to that of a totally free theory. So 
although the eigenvalues may not accumulate fast enough near A = to produce a chiral condensate, 
the accumulation rate may still be very different from that of a theory with no bare interactions at 
all. This raises the question whether also such spectral densities display universal features that are 
computable from random matrix theory. As mentioned, it is even quite likely that QCD possesses a 
large "conformal window" in a range of flavors, where Nj < 16 but still greater than some critical 
value A^r. In such a phase, which is not expected to support the spontaneous breaking of chiral 
symmetry, the spectrum of the Dirac operator will definitely not be that of a totally free theory 
either. Finally, although the scenario is quite different in three dimensions (because the gauge theory 
is superrenormalizable there, with a trivial beta function), the outcome of adding more fermionic 
flavors to the theory is believed to be quite similar. More flavors imply more screening of charges, and 
hence reduced forces between the fermions. Eventually the forces should become too small to cause 
chiral symmetry breaking. So even if flavor symmetries are spontaneously broken in three dimensions 
for a low number of fermionic flavors, one does not expect this to persist as Nf is increased. 

Yet another example may be general SU{Nc) gauge theories coupled to both fermions and Higgs fields. 
In the Higgs phases of such theories we again do not expect chiral symmetry breaking. In the absence 
of Yukawa couplings the Dirac operators of such theories are identical to those without Higgs fields, 
and the arguments of ref. ||l| would naturally lead one to associate analogous matrix models with 
vanishing spectral density at the origin to such phases. In fact, if one mentally integrates out the 
Higgs fields, the resulting theories are formally of the same class as the theories without Higgs fields. 
Their only difference is that they belong to different phases. A totally different class of gauge theories 
in four dimensions which do not cause spontaneous chiral symmetry breaking are those of Abelian 
gauge groups, such as QED. Also for these much simpler theories one can ask whether their Dirac 
operator spectrum near the origin contain universal functions. 

All of these physical considerations have led us to analyze the fate of the universal eigenvalue correlators 
(and in particular their microscopic spectral densities) when the the macroscopic spectral densities 
at the origin reach zero, i.e. when p(0) = but without a finite gap in eigenvalues around A = 0. 

^ The number Nf = 17 at which chiral symmetry breaking no longer occurs is probably too conservative. Lattice 
gauge theory simulations indicate that it may occur much before pO[. 



Prom a pure matrix model point of view, these questions are clearly posed, and may be considered 
in their own right. For example the proof for the macroscopic universality of density-correlators with 
such a multi-band support for the eigenvalues has been recently given for the hermitian and complex 
matrix model |]l^]. Their possible physical relevance in terms of applications to the Dirac operator 
spectrum near the origin can necessarily only be phrased in terms of conjectures. They are, however, 
highly non-trivial statements, and they are directly testable in, for example, lattice gauge theory 
simulations .0 

To require that the (macroscopic) spectral density precisely just vanishes at the origin puts a constraint 
on the matrix model in question. This implies the fine tuning of at least one parameter (coupling) in 
the matrix model potential V{M). From this point of view it is not surprising that if we associate the 
usual scaling (with no such constraint) with normal critical behavior, the just- vanishing of the spectral 
density at the origin corresponds to, in general, multicritical behavior. There are, for similar reasons, 
tight constraints on the precise manner in which the macroscopic spectral density can approach zero 
at the origin, in general labeled by an integer m, which counts the order of multicriticality. To achieve 
such multicriticality of mth order, m couplings must be tuned. This integer m in turn describes the 
rate at which eigenvalues accumulate near the origin, a rate which hence cannot be totally arbitrary. 

With these preliminary remarks to serve as motivation, we now turn to the details of the actual 
matrix model calculations. Our analysis will be based on the orthogonal polynomial technique, which 
is known to greatly simplify many calculations in the unitary and chiral unitary matrix ensembles. 
As in ref. 0, we shall turn the three-step recursion relations for these orthogonal polynomials into a 
differential equation in the large- iV limit. Prom this we can derive the limiting forms of the orthogonal 
polynomials in the appropriate large-A^ scaling regions. A recent paper by Kanzieper and Freilikher 
p^ shows very elegantly how the derivation of this differential equation can be simplified. Their 
approach is in fact more powerful, and leads more easily to generalizations. We shall essentially follow 



their approach here. As derived in ref. |13|, the crucial property of the differential equation is that 
it does not explicitly depend on the chosen matrix model potential V{M) once the large- A^ limit is 
taken, but only on the associated macroscopic spectral density p{X). As we shall see in section 2, 
this simplifying feature unfortunately does not survive in the particular microscopic large-A^ limit 
that is required to probe the multicritical domain near A = 0. The resulting differential equation will 
therefore not manifestly depend only on the macroscopic spectral density p{X), but will also contain 
contributions from what are naively subdominant terms in the 1/A^ expansion. Universality of the 
orthogonal polynomials, and hence universality of the microscopic spectral density and the spectral 
correlators, is hence not as obvious as in the usual case. Finding those naively subleading corrections 
that turn out to contribute in the microscopic multicritical limit entails solving an auxiliary differential 
equation, known in a different context as the string equation.^ For the case at hand it turns out to 
be a generalized Painleve II equation, which we show how to solve accurately by numerical methods. 
Both because the solution to the generalized Painleve II equation cannot be found in closed analytical 
form, and because the underlying differential equation for the orthogonal polynomials themselves 
cannot be recast in a particularly simple form, it is not possible to give a closed analytical expression 
for the microscopic spectral correlators in this multicritical case. The solution can be found to any 
required accuracy by numerical means. We do however identify an intermediate "mesoscopic" range at 
which the spectral correlators can be computed analytically, expressed by some combination of Bessel 



■^For a beautiful comparison between recent lattice gauge theory data and the universal predictions from matrix 
models in the usual case of chiral symmetry breaking, see ref. |l3] . 

■^Only after properly solving this additional differential equation, and inserting its solution, does the differential 
equation for the orthogonal polynomials become manifestly universal. 



functions. In section 3 we derive, for completeness, the analogous results for the chiral-unitary matrix 
ensembles. They may be relevant for Dirac operators of fermions in the fundamental representations 
of gauge groups SU{Nc'>3) in even space-time dimensions. These results follow as simple corollaries 
from the unitary ensemble. Section 4 contains our conclusions. 

2 Multicritical Unitary Ensembles 

We begin with the unitary matrix ensemble described by the partition function 

Z = /dMdet^^M e-^^^^(*^) 

N 



r n (dA.IA.p-e-^^^^^)) det.,A;.-i ' , (2.1) 



with some even potential 



V{M) = J2 ^M'^^ . (2.2) 

The matrices M are N x N and hermitian, with real eigenvalues A,. The integration is over their 
associated Haar measure dM. This matrix model is, for generic potentials V{M) giving rise to a single 
cut containing the origin of eigenvalues A = 0, conjectured to be related to the Dirac operator spectrum 
of SU{Nc'>3) gauge theories coupled to an even number of flavors Nf = 2a. For that application a is 
thus restricted to be an integer. Using the formalism of Kanzieper and Freilikher [14| (see also refs. 



|15| , |lq| ) this restriction can be lifted without difficulty, and we shall therefore in greater generality 
consider any real a > —1/2. We shall largely employ the notation of ref. [0] in what follows below. 

Using the method of orthogonal polynomials to analyze the above partition function, one has to find 
polynomials -Pn(A) orthogonal with respect to the weight function w{X) = |Ap" exp[— A^1/(A)]. We 
choose them orthonormalized: 



dX w{X)PmiX)Pn{X) = 5mn ■ (2.3) 

D 

It is actually more convenient to work directly with the "wave functions" 

V'„^(A) ^ |Are"(T)mp^(A) , (2.4) 

because these are the objects that enter directly into the determination of the spectral kernel (two- 
point function): 

N-l 

Km{X,Ii) = ^ tpn{X)tlJn{fJ') 

n=0 

V'Af(/^)V'Ar-l(A) - V'Ar(A)V'Af-l(M) /„ PN 

= CAT . (2.5) 

/i — A 

Here the coefficients cn are those of the three-step recursion relation for the orthogonal polynomials: 

AP„_i(A) = c„P„(A) + c„_iP„_2(A) . (2.6) 

They can be determined from the following identity^ (which in the context of 2-d quantum gravity is 



called the string equation) by partial integration, using eq. ( |2.6[ ), 

POO 

n = cn dXw{X)P^\X)Pn-i{X) . (2.7) 



*By switching to the monic normahzation P„ (A) = A" + . . . one can easily derive that P„' (A) = —P„-i{\) + 0(A" 



The wave functions V'm(A) are clearly orthonormalized on the usual (flat) Lebesgue measure. 

In order to derive a differential equation for the wave functions it is convenient to express the differential 



operator acting on the orthogonal polynomials in the following way [14|: 



P„'(A) = ^„(A)P„_i(A) - 5„(A)P„(A) , (2i 



where the functions An and i?„ have to be determined. Following Kanzieper and Preilikher |L4| , they 
are given by 

^n(A) = ^^c„|_'2d^^[F('„)(^)-y('„)(A)]P.(/.)^ (2.9) 

BnW = iVc„ rdM^[^(;)(/^)-^H(A)]Pn(/i)Pn-l(/^) 

= i?-s(A) + (l-(-l)«)^, (2.10) 

where the a-dependent part of the measure has been absorbed into the potential^: 

y(,)(A) ^ F(A)--ln|A|. (2.11) 

In contrast to reference [jl^ we have chosen the usual matrix model convention eq. (^]^) which makes 
the dependence on A^ explicit. The evaluation of different powers in N will become very important 
in the multicritical case. Using the results eqs. (|2.9|) and (|2.10| ) as well as the recursion relation eq. 
( ^ ) the following identity is exactly satisfied Q: 



Bn{\) + Bn-i{\) + NVl^^{\) = -^An^i{\) . (2.12) 

C-n—1 



Differentiating eq. ( p.SD for Vra(A) and using relations ( |2.6D and ( 2.12 ) finally leads to the following 



exact differential equation for the wave functions V'n(A), valid for any n |L^]: 

<(A)-F„(AX(A) + G„(A)V^„(A) = 0, (2.13) 

where 

F„(A) . fl (2.14) 

G„(A) s -2^/l„(A)/l„_i(A)-fB„(A) + :^V,'„,(A) 

+ (S„(A) + yF(',)(A))' - ^ (i?„(A) + ^Vl^^{\)) . (2.15) 

The rest of our paper will consist in analyzing the way in which suitable microscopic large- A'^ limits 
affect the above differential equation. Performing the usual large-A^ limit the recursion coefficients 
Cn and the functions A„(A) and B^^^{X) are assumed to approach smooth functions. Under these 

^In An{X) the term proportional to a vanishes due to the potential V(c,)(A) being even, whereas for B„(A) it is given 
explicitly, following again pM. 



conditions, with a generic potential V{X), the function An{X) is directly related to the macroscopic 
spectral density p{X) |jl3|, ^]: 



P(A) 



hm — A;v(A)v^l-(A/a)^ 
N^oo iVvr ' 



(2.16) 



where a = limAr^oo cat is the endpoint of the spectrum. Consequently, because of eq. ( 2.12| ), the 
asymptotic differential equation ( p. 13 ) for tpj\[{X) should only depend on the macroscopic spectral 
density. But precisely at multicriticality, and in the appropriate microscopical scaling limit, the 
simple relationship ( p.l6| ) unfortunately breaks down. This is due to the contribution of terms that 
are naively subleading in 1/A^. Because of this, the appropriate differential equation for the wave 
functions ^l^]s[{X) in the multicritical microscopic limit will not just depend on the spectral density 
p{X). To see this, we first make the definition of multicriticality precise. 



2.1 Multicritical potentials 

Multicriticality can be reached at both the endpoint of the spectrum, and at the origin. Because of 
the physical motivations mentioned in the introduction, we here focus exclusively on multicriticality 
at the origin of the spectrum. The case of p{0) ^ will thus play, in this context, the role of being 
uncritical {m = 0). As p{0) -^ 0, the one-cut solution around the origin will turn into a two-cut 
solution, thereby defining the first criticality {m = l). As is well known, this happens in the simplest 
case ll^, ll^, pO] when the potential V{X) has a double-well shape, as shown in fig. 1. 
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Figure 1: The critical m=l potential (left) and the associated spectral density (right). 
A characterization of potentials leading to multicriticality at the origin has been given by Crnkovic 



and Moore |21|. In general, for polynomial potentials of degree 2k, the macroscopic spectral density 
is, for a symmetric one-cut solution, of the form 



p{X) = Q2k-2{X)^l-{X/af 



(2.17) 



where Qk{X) is a polynomial of degree k. Generically, Qk{0) 7^ 0, in which case p{X) ~ const, near 
A = 0. To force p{0) -^ thus entails tuning the potential V{X) so that the polynomial Q{X) vanishes 
at the origin. Because Q{X) is a polynomial of even degree, this can be achieved by powers of an even 
integer, 2m, where m labels the order of multicriticality: 



P(A) 



const. A^™ + 



as A — > 



(2.18) 



Polynomial potentials Hn(A) that lead to such behavior are thus at least of degree 2m+2. A particular 
class of such potentials for which the behavior of the macroscopic spectral density is as above can be 
compactly described by |21]j 



y^(A) = k{m)\ 



2m+l 



1 



1/2 



(2.19) 



where endpoint of the support a has been set to unity and the "+" subscript indicates that only 
the polynomial part in an expansion around A = oo is kept. The constant in front, chosen so as to 
normalize the macroscopic spectral density by 



d\p{\) 



1 



is given by 



k{m) 



,2m+i(rn±mrn-i)\ 



(2m -1)! 
The spectral density at the multicritical point is then given by 

1 



/Om(A) 



2tt 



k{m)X^"'Vl-X^ . 



(2.20) 



(2.21) 



(2.22) 



Conforming with common usage, we shall denote such multicritical potentials as "minimal". An 
infinity of other potentials can be constructed by adding higher powers of A the polynomial Qfc(A), 
which lead to precisely the same multicritical behavior ( 2.18| ). For this reason one can speak about 
an mth multicritical universality class. One sees immediately that the a-dependent term in our 
integration measure has no influence on the classification of these multicritical potentials, since this 
term is suppressed by one power of 1/N in the planar limit (see eq. ( p.llD ). 



2.2 The breakdown of the usual planar limit 

As it has been mentioned earlier in this section an asymptotic differential equation for 'il^]\f{X) at 
multicriticality is no longer easy to obtain due to a breakdown of the relation ( 2.16| ). In order to show 
this let us briefly review how the result is obtained in the noncritical case, where the usual planar 
limit produces all relevant terms, and point out what breaks down at multicriticality. 



In the usual planar limit the following behavior is assumed for the recursion coefficients: 

Cn±l ^ Cn + 0{1/N) , 



(2.23) 



which is directly related to the fact that we are here considering a one-cut phase. Moreover, an implicit 
assumption in ref. |14| is that the following smooth limits are satisfied: 



A„±i(A) ^ AniX) + Oil/N) , 
Bl'l,{X) - B^:^{X) + 0{1/N). 

The requirement eq. ( |2.23| ) together with the leading term of the relation 



2CN 



a + ... 



(2.24) 
(2.25) 



(2.26) 



Crnkovic and Moore have considered a spectral density with two symmetric cuts merging. This leads to the the 
same multicritical potentials Vm(A) as a symmetric one-cut density developing zeros at the origin. 



between the coefficient cn and tlie endpoint of the support leads to the asymptotic relation ( 2.16 ) 
between An{X) and the macroscopic spectral density p{X). Inserting the limit ( |2.25 ) into eq. (2.12) 
leads asymptotically to 

Bn{\) + ^VU^) - ^^n{\) - (-1)^^ . (2.27) 

Consequently GnW can be expressed asymptotically by Aj\f{X) (and via eq. ( |2.16| ), if true, by p{X)) 
only[|. 

G^{X) ^ ^^(1 - ^) + -^ + (-1) --^-^ . (2.28) 

Looking at the noncritical case at the origin the terms proportional to -j^^-jx) ^^^ subleading in the 
microscopic limit at the origin so the last term in Gat (A) and the -F/v(A)-term can be disregarded. 
This leads to the asymptotic differential equation 

i^Nix) + ( ^^2^ + ^-^2 1 V'iv(x) = , (2.29) 

where the appropriate scaling x = NX has been used and Va^(^ ~ ;^) ^^ regarded a function of x. 
It has been derived and solved for nonnegative integers in 0, and generalized to real a > —1/2 in 
jl4{ |. Apart from the independence of the two derivations, we have also checked the resulting solutions 
(certain Bessel functions, see ref. Q) by numerically determining ijjn{X) directly from the recursion 
relations themselves, for very large values of A^. 



Moreover there is an analytic consistency check on the differential equation ( 2.29| ). It is straightforward 



to verify (see Appendix A) that the following exact symmetry holds for the wave functions: 

4r^)(A) = sign(A)4:Vi(A) ■ (2.30) 

The superscript indicates the corresponding value of a in the measure. We will refer to this symmetry 
as a-symmetry. Note that the macroscopic spectral density is independent of a, since the a-dependence 
can be considered as a subleading term in the potential. Due to relation (2.16) and the limit ( |2.24 ) 



the differential equation (2.2S|) is therefore manifestly invariant under the a-symmetry. Indeed, the 



solutions of ref. [^ trivially satisfy the relation ( 2.30| ). 



Let us now turn to the multicritical case. The appropriate rescaling of eigenvalues that is needed to 
obtain a non-trivial microscopic limit will now be of the form 



1 



Af^^^^A (2.31) 



at the mth multicritical point (m G N). We first suppose that the usual planar limit eqs. ( p. 23 ) - ( |2.25D 



still holds. One might object that the recursion coefficients now obey a period-two ansatz |19| instead 
of eq. ( ^.231) . But in the planar limit, and just at the critical point n = N, the two branches merge 
to the one-cut value cn = a/2, as we are right on top of the one-cut-two-cut transition. Therefore 
the argument in |14| relating A]\j{X) and p{X) (eq. (|2.16D ) should still be expected to hold. Next, we 
have seen from the previous subsection that at the mth multicritical point p{X) behaves like A'^™' at 



A = 0. Consequently the last term in eq. ( 2.28| ) for Gat (A) becomes {—)2ma/}? which is of the same 
(leading) order as the second term. 

^The last term has been omitted in eq. (34) of ref. |l4| because the authors did not consider the microscopic limit at 
the origin. 



So from the usual planar limit we are left with the full Gat (A) eq. ( |2.28| ) in the asymptotic differential 
equation for ■07v(A). However, this differential equation is no longer invariant under the a-symmetry. 
It therefore cannot be the correct equation. In order to confirm this, we have compared the prediction 
from this naive differential equation to the polynomials constructed by Gramm-Schmidt orthonormal- 
ization. The mismatch is shown in fig. 2. 




>t%+i [x] 




Figure 2: The microscopic wave functions for the m — 1 critical potential ( ^.1S| ) are plotted for A^ = 32 and 48 
(a — 0) (real lines). They converge to limiting functions which differ from the solution of the naive differential 
equation (^.28) (dotted lines) for small x. 



The only explanation for this failure is that the simple relationship ( |2.I6| ) breaks down or, in other 
words, that the a-dependent corrections to the recursion coefficients that were subleading before have 
to be taken into account. Remarkably, this can also be established directly by analyzing the "string 
equation" for the coefficients c„, as we will discuss in section 2.4. One finds that corrections to the 
c„'s are precisely of the order N~'^ for 1^= l/{2m + 1) at mth multicriticality, instead of N^^ (as for 
the case m = 0). It is due to this rather devious compensation that the naive differential equation for 
the wave functions is incorrect. 



Despite many efforts we have not been able to derive a general expression for the asymptotic differ- 
ential equation containing all relevant contributions from the CnS. The analytical formulas increase 
dramatically in complexity with increasing m. We will therefore restrict ourselves to a detailed ana- 
lytical treatment of the m = l and m = 2 multicriticalities. Many aspects of the formalism do however 
readily generalize to arbitrary m, and we shall keep the discussion in as general form as possible. In 
Appendix B we write down a detailed conjecture for the differential equation for the wave functions 
at any multicriticality. 



2.3 An EXAMPLE: THE CASE m = 1 

In order to treat all contributions to the differential equation at multicriticality we will restrict ourselves 
first to the simplest example, the case m = l. This corresponds to a minimal potential of the simple 
Mexican-hat form 

Vi{X) = ^52 A^ + -54 A^ with 52 = -8 , 54 = 16 . (2.32) 



The corresponding eigenvalue density reads| 



zvr 



(2.33) 



We are now able to give exact expressions for the functions An{X) and Bn{)^) only in terms of the 
recursion coefficients Cn- This will allow us to keep track where the correction to the Cn's will enter. 
We will then be able to derive an asymptotic differential equation at tti = 1 multicriticality. 

Making use of the completeness of the orthogonal polynomials we can explicitly express -P„'(A) by 
lower order polynomials, as we have chosen a particular potential: 



/oo 
dlXw{p.)NVl^^{^l)Pn{lj)Pk{lj) 

= NCn (52 + 54(4+1 +ci + ci_S) Pn-liX) + iVc„C„_iC„_294i^n-3(A) - (1 - (-l)")yPn(A) 



(2.34) 



where the recursion relation (2.6) has been used for the potential eq. ( p. 32 ). Looking back at the 
definition ( 2.34| ) An[X) must be an even function whereas B^^{X) is odd. A simple way of finding them 
is to assume that they are polynomial^ use the recursion relation and then comparing coefficients. In 
this way one finds 



An{X) 
BniX) 



Ncn [92 + 54(0^+1 + c^) + g^X- 
Nclg^X + (l-(-l)")^. 



(2.35) 
(2.36) 



One can also derive these expressions directly, using the the definitions ( |2.9D and ( 2.1C| ) . If we would 



now set c2_^i = c^ : 



1/4 and use the critical couplings eq. ( 2.32 ) the combination 52 + 5'4(c^+i + c^) 



would vanish and we would be back at the relation ( ^.16 ), which no longer holds at multicriticality. 
From our example we are now in a position to determine the order and the explicit value of the 
corrections. Inserting the appropriate microscopic scaling limit at multicriticality, eq. (2.31), the 



whole differential equation ( 2.13| ) has to be rescaled by A^ 2m+i ^ which is N ^'^ in our case. In order 



to get contributing corrections in the differential equation, the following terms have to be of order one: 



V = m (52 + 2544) . 



(2.37) 



In other words, for m = l the terms in the parenthesis must vanish like N~'^i^ and A~^'^, respectively. 
This is precisely what we get from eqs. ( 2.44 ) and ( ^.51 ) below. Expressed in terms of the above 
quantities, we get the following exact differential equation for the potential eq. ( 2.32| ) (from eqs. 
(HI), O) and (|1D): 



^'^{x) - N-^'^Fn{x)^'m{x) + N-^I^Gn{x)iI:n{x) 







(2.38) 



®Here a = 1 due to our conventions from subsect. |2.l| . 

^Remember that in the usual planar limit An{\) becomes proportional to the polynomial Qk{^) in front of the square 
root of p(A). 
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Here 

N-^/^Fn{x) 



n+ + 54X2 



N-^I^Gn{x) = c%u+u^ + ((-l)^a + ly+ (^^(^+ + "-)^4 - \v^ + iV-i/3((-l)^a + ^)54) x 



+ ( c%al - '-N-y^vg,) X- - \N--I^glx^ + i-^r^^.a-vg^x^-N-^^glx^ 



2/4 u+ + ^'40;^ 

+ ^ ^ " . (2.39) 



X 



Due to eq. ( |2.37| ) it is now clear that the terms being explicitly proportional to negative powers of A^ 
will vanish in the microscopic large- A'^ limit. Moreover, the dominant term ~ x'^ in G]\j{x) vanishes 
up to 0(Ar-i/3). 

v'^ = Ac%g^{u+ + u^) . (2.40) 

This follows directly from the string equation ( p.Tj ) at n = A^, which in our example reads 

^ + (1_(_1)-)^ = c2 (52+ 54(4+1 + c2+c2_i)) , (2.41) 

and from inserting the critical values for the coupling constants. Finally the correct asymptotic 



differential equation for m = 1 criticality is eq. ( 2.3§| ) with 



^T ■?n^ / N 9 f, -.sN lA u+f + (— 1) 2054 994 i—l) a — a 



2 I '-JVi'4-^ I 9 

u-^- + (74X^ x^ 

(2.42) 

where one of the quantities u-^,u- or v can still be eliminated by means of eq. ( 2.40 ). After having 
determined the functional form of the differential equation, which is now no longer a Bessel equation, 
the only point left to be done is the analytical calculation of the corrections to, say, c^_|_i and c^. This 
will be the topic of the next subsection. 

Let us add a final remark on the a-symmetry of this new asymptotic differential equation. The 
function Fn{X) can be seen to be invariant for finite n from the definition eq. ( 2.14| ) since one has 



4r^^(A) _ 4ii(A) .,43^ 



This follows directly from eq. ( 2.30| ). An immediate consequence of this identity is the invariance of 



M+ under the a-symmetry. Still, apart from the last two terms, G]\f{x) is not obviously invariant. In 
order to prove its invariance under the a-symmetry one would have to know precisely the behavior of 
U- and V as functions of a. As we shall see below, this is a highly non-trivial problem, related to the 
solution, for arbitrary a, of a certain non-linear differential equation. 

2.4 Painleve II 

The apparently subdominant terms that enter the differential equation for the wave functions ip{X) in 
the microscopic limit near multicriticality arise from corrections to the coefficients c„ as A^^oo. It 



is well-known |19, 20| that in the case at hand the coefficients Cn do not approach one single smooth 
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function in the large-A/" limit, but rather spht into two functions depending on whether the index is 
odd or even. In ah generahty pI| ] P^ 

cln = C + N-^fmiz) + N-^'^gmiz) + N-^'' h^{z) + ... 
cln+1 = C + N-'^frniz-N-'^) + N-^''g^{z-N-'') + N'^'^kniz - N~'') + ...(2.44) 

for 2n ^ N, where the scahng variable z is defined by 
with the exponent 



z = iv^™^(l-_l (2.45) 



1 



(2.46) 
2m + 1 ^ ^ 

at TTith multicriticality. Here fm{z),gm{z), hm{z) and the corresponding barred quantities are assumed 
to be smooth functions. 

Consider first the case m = l. At the critical point, which is at -j^ = 1, the two branches for the 
even and odd recursion coefficients meet in the planar limit, which implies C = C. Inserting eqs. 
(p. 44) into the string equation (|2.41|) for even and odd n, and expanding consistently in N~'^ up to 



third order, yields the following set of equations: 



C = C = \, (2.47) 

which is nothing else than the statement 2cn = a + . . . . Furthermore we get 

fi{z) = -h{z) (2.48) 

gi{z) + -gi{z) = 2h{zf-^z (2.49) 

= /i(z)(l6/i(z)2-z)-2/('(z) + |, (2.50) 

where the latter is a Painleve II equation |2^] for a = 0. The dependence on hi{z) and hi{z) has been 
eliminated using the equations from even and odd n. Due to the following observation [^] 

cL+, + cL+g+i = 2C + N-^m-lYC{z)+gm{z) + -gm{z)) + ... , (2.51) 

where 2n + q ^ N, we just have to determine the combination gi{z) + gi{z) and the function fi{z) in 
order to get the quantities u± and v from the last subsection. So we just have to solve the generalized 
Painleve II equation ( ^.50 ). 



Consider first the case a = 0. The solution to eq. (2.50) can then be shown to be unique |20, ^ once 



the following physical boundary conditions are imposed. First, in the usual planar limit z -^ oo, we 
takep*] fi{z) = — -y/z/4. This follows directly from solving the string equation ( |2.4lD in the large- A'^ 
limit. Next, we demand that there are no singularities on the real z-axis. (This is equivalent to 
enforcing fi{z) = at z = — oo). Because the boundary conditions in this manner are given at the 
opposite ends of the interval (here the real line), and because the solution is unstable towards the 
development of singularities once slightly improper boundary conditions are imposed, it is not entirely 
trivial to compute the correct solution numerically. We have found that the discretized, iterative, 

^°Here we have used the most general ansatz for the coefficients c„. 

^^With fi{z) being a solution also —fi{z) is a solution. We fix the sign by noting that for any finite A'^ the branch of 
the coefficients c^ with odd n is always larger than the one with even n |23]. 
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procedure proposed in ref. ||2J] for the odd-order generalization of Painleve I equation works well 
here too. In fig. 3 we show a plot of the solution to eq. ( 2.5C| ) found in this way. As a check, we 
display on the same figure a numerical determination of fi{z) directly from the recursion coefficients 
c„, evaluated to high order. For detailed definition of the parameters used to obtain this numerical 
solution, we refer to ref. [p^j . 



-f [z] 



-4 



uuih:!":"" 

-2 



0.5 



0.3 



0.2 



a::i 



Figure 3: Numerical solution to Painleve II equation f{z){16p{z) — z) — 2/"(z) = (lower dots). We have used 
the parameters: cutoff= 5, mesh= 1/10, calibration parameter h — 1/300, number of iterations v = 400 ~ 500, 



and extrapolated the result to i^ = oo. Explicit values of /„ = iVa (2c^ 



^n+l 



/4 for the critical 



potential (2.32) with a = Q, N = 1024 are also plotted (upper dots). 



Inserting (72 = ~8 and 5(4 = 16 for the minimal m = \ multicritical potential, we get for A^ evenP| 



n± = 16[2/i(0)2±/{(0)] 
V = 32/1 (0) . 



(2.52) 



Note that the leading contributions (proportional to N^ and As, respectively) to the above quantities 
vanish at the critical point. The first corrections, from eq. ( |2.51| ), are precisely of the required power 
to render both u± and v finite in the A ^00 limit. Yet higher order corrections can indeed be ignored 



here. If we substitute the solution to the Painleve II equation, the full differential equation ( 2.38| ) for 
the orthogonal polynomials in the microscopic limit is thus known in this case. Numerically we find 
from the above solution u+ ~ —1.105, u- ~ 2.211 and v ~ —4.207. In fig. 4 we plot the solution to the 



equation (2.38) obtained in this way. On the same figure we also show the wave function for ipi^{x) 
constructed directly by Gramm-Schmidt orthonormalization. The two functions are seen to agree to 
very high precision. This gives us complete confidence in our procedure. 



"For TV odd one has to replace ± -^ =p on the r.h.s. and v 
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Figure 4: Numerical solutions to the linear differential equation (2.3S) using data from Painleve II (bold 
lines) . Their vertical scales are chosen to optimize the matching with the microscopic wave functions {N = 48) 
constructed from the critical potential with a = by Gramm-Schmidt orthonormalization (real lines) . 

The microscopic correlation functions are constructed out of these limiting functions for the wave func- 
tions. In fig. 5 we exhibit a microscopic spectral density ps{x) ~ ('i/'Ar(x)V'7v-i(^) ~ tpN-i{x)i^'^{x)): 




0.25 0.5 0.75 



1.25 1.5 1.75 



Figure 5: Microscopic spectral density for the m — 1 criticality, divided by x^ in order to take out the trivial 
x^-growth which is a consequence of the matching-on to the macroscopic spectral density. The normalization 
is determined by requiring \imx ^ oo p's{x) = p"(0) = 16/7r (dotted line). 

Next, consider e.g., the case a = l. Appropriate boundary conditions for the generalized Painleve II 
equation are easily established. In the planar limit, here reached at z — > ±00, the a-term does not 
contribute, as it is suppressed by 1/A^ compared with the rest of the potential V'(A). The boundary 
conditions are therefore precisely the same as for a = 0: 



fi{z) — > as z ^ —00 , 

fi{z) — > —y/z/A as z — > 00 . 

For numerical purposes more refined boundary conditions are convenient: 



h{z) 



al{^z) 
1 / z 



as z 



-co 



.3y{z 

with y{z) 



+ y{z} + 



iVs r 



-a + \ a 



\3y{z) 

^3/27 



viz) 



as z 



00 



(2.53) 
(2.54) 



(2.55) 



(2.56) 
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This last condition is found from solving the algebraic 3rd order equation obtained from eq. ( |2.50D , 
self-consistently setting fi{z) ~ as z — > oo, and choosing the root which reduces to fi{z) = —^/z/4: 
in the limit z—>-oo. In fact, because the solution (|2.54D , extended to the negative z-axis by 



^'^'^ = vkz)^\y^'^ 



for 2: < 



(2.57) 



has /('(-z) ~ everywhere, and automatically satisfies boundary condition ( |2.53| ), it is almost exact. 
Numerically it is therefore an extremely efficient starting point for finding the full solution to eq. 
( p. 50 ). It becomes exact in the limit a ^00 (see also below). 



In fig. 6 we plot the solution of the differential equation ( p. 38 ) for a = 1, and compare it with the 
solution explicitly solving the recursions relations, using very for large values of N . We again find 
perfect agreement. 
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Figure 6: Numerical solution to Painleve II equation /(2:)(16/^(z) — z) — 2f"{z) + 1/2 = (lower dots). We 
have used the same parameters as in fig. 3. Explicit values of /„ for the critical potential with a = 1, A^ = 512 
are also plotted (upper dots). 



This is also a good point to check that the a-symmetry ( |2.30| ) is satisfied once we take the solution 
to the generalized Painleve II equation into account. It follows from eqs. ( p.30| ) and (2.50) that the 
required relation is (indicating the a-dependence explicitly by a subscript): 

/{(,) (z) + 2 A(,) {zf = - /{(,_y (z) + 2 A(,_i) {zf . (2.58) 

This is obviously a highly non-trivial identity which the appropriate solution to the generalized 
Painleve II equation ( 2.50D must satisfy. In the limit a — > 00 the simple algebraic solution (2.56) 



approaches, for fixed z, the full solution to the generalized Painleve II equation ( |2.50p . It is straight- 
forward to verify that the algebraic solution (|2.56| ) indeed satisfies the identity (|2.5^ ), order by order 
in a large-a expansion. Moreover, by performing two Miura transformations we can actually prove 
that the relation ( |2.58| ) is exact. The two transformations are 

M{z) = f[{z) + 2h{zf (2.59) 

M{z) = -f[{z) + 2fi{zf and shift a^a-1. (2.60) 

Under these two different transformations the generalized Painleve II equation ( p. 50 ) turns into the 
same equation: 



a(l -a) + 2z'^M{z) - 32zM{zf + l28M{zf - AM'{z) + 16 [M'{z)\ '^ + 4zM"{ 



■32M{z)M"{z) = 0. 
(2.61) 
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The two functions ( 2.591) and ( |2.60| ) thus satisfy the same differential equation. The two boundary 



conditions are of course inherited from those of fi{z), and one can readily check from eqs. ( 2.53| )-( |2.54 ) 



that the two functions ( p.59| ) and ( 2.60| ) also satisfy the same boundary conditions at z^itoo. The 



two functions must therefore be identical, and the relation ( 2.58| ) has been proven. 



Note that, to the order at which we are working, the identity ( 2.581) implies the following relation 



among the recursion coefficients (indicating the a-dependence explicitly by a superscript): 

(a)2 ia)2 _ (a-l)2 (a-l)2 ,„ fi2^ 

'^2n "T ^2n+l — ^2n+l ">" ''2n+2 • \^.UZ,) 



In this form one easily proves that the differential equation for the wave functions ( 2.38| ) respects the 
a-symmetry ( 2.30| ). 



2.5 Universality 

Up to now we have only been investigating the most simple example, m = 1 multicriticality with the 
minimal potential. In order to prove universality, one has to show that all perturbations to the m = \ 
minimal potential that remain within the m = l class lead to the same universal differential equation 
( ^.38 ). It should be parameterized by only one universal constant, /o"™''(0), in direct analogy with the 



noncritical case m = in [Q]. Since there is no longer a simple relationship as eq. ( p. 16 ) at hand, this 



is a rather involved task. Still, we have strong evidence that universality holds precisely in the sense 
just described. In order to illustrate this, we will give as an example a nonminimal m = l multicritical 
potential. We will also briefly outline the situation for m>2 below. 



When adding a sextic term to the minimal potential eq. ( 2.32| ), 



Fi(A) = ^g^X^ + ^ff4A^ + ^<76A^ , (2.63) 

the conditions for m = 1 multicriticality and the corresponding eigenvalue density read 

= 52 + -54 + -56 , 

1 = ^(54+56), (2.64) 

16 
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P(A) = ^(^(54+256)A' + 56A^j V1-A2. (2.65) 

The outcome of this modification is that the calculations of the previous 2 subsections exactly go 
through when replacing 

54 ^ 5* = 54 + 2^6 , (2.66) 

or, in other words, it depends only on one single parameter ~ p"{0)- 

Without going through the details, where one can use the results for m = 2 from appendix^, we define 
u± and V as in eq. ( |B.5D , but with a scaling in powers of A^ as in eq. ( 2.37| ). The final result reads 



u± = 5*(2/'(0)±/'(0)), 
V = 2gJ{0) , (2.67) 
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where the generahzed Painleve II equation is now given by 



= gJizf - zf{z) - g^-fiz) + | . (2.68) 



Using eq. (|2.67|) and g^ instead of (74 in eq. ( 2.42| ) for Gn and eq. (|2.3g|) for F/v provides the 



differential equation ( |2.38| ) for the m = 1 nonminimal potential eq. ( |2.63| ), where (7* = p"{0)/Tr. Eqs. 
( ^.67 ) and (|2.68| ) and thus the differential equation ( 2.3§| ) with (|2.42| ) will hold also for all higher order 



perturbations of the minimal potential. 

The essential ingredient in obtaining universality is thus the fact that we still have the freedom of 
modifying the potential while remaining in the same universality class. The only quantity that enters 
in the differential equation is, in the m = 1 case, the second derivative of the macroscopic spectral 
density at the origin, p"{0). All other details of the specific m = l potential have disappeared. By a 
simple rescaling of p"{0), all m = l critical potentials should thus lead to the same universal functions. 

The way this should generalize to the higher-m case is as follows. An alternative way of characterizing 
multicritical potentials is by the scaling behavior of the recursion coefficients c^ near the "critical 
point", where the two branches of the period-two solution merge. Here, in the usual microscopic 
limit H], 

/-(^) = - (22n.+i(^ + i))iA>n ' for z - 00 (2.69) 

which is the generalization to higher m of the behavior lim^-^oo /i(-z) ~ ~\f^l^ we found explicitly 
in the m = l case. This scaling of the recursion coefficients near the point where they merge, cf. eqs. 
( p^ and ( ^Wj ), 
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4-C ~ ± const. (^1-^j (2.70) 

is universal in the sense that its power-law scaling depends only on the class m of multicriticality (the 
positive solution corresponds to odd n, while the negative one corresponds to even n). The coefficients 
C = C' of eq. ( [2.44]) will change with the chosen 771th multicritical potential, and so will the particular 



normalization of (2.69|), but the scaling (|2.70|) will remain. The non-universal coefficients C = C drop 



out of the differential equation for the orthogonal polynomials simply by the requirement of mth 
order multicriticality. The next term of the expansion ( p. 51 ) is what will enter in the differential 



equation. Our statements about universality above translate into saying that for all m this particular 
combination of solutions to the generalized Painleve II equations for fmiz) (and higher functions of 
the expansion ( 2.44| )) will be proportional to p^'^^'{0). We have not proven this. 



2.6 A Double-Scaling Limit 

As we have shown above, apparently non-leading terms from the coefficients c„ turn out to contribute 
in the multicritical microscopic limit. At the core of this phenomenon stands the breakdown in this 
limit of the simple relation (p. 16) between the function An (which enters directly into the differential 



equation for the polynomials), and the macroscopic spectral density p{X). The way this relation breaks 
down at multicriticality makes it clear that the phenomenon is more general, and directly related to 
the very definition of multicriticality. So far we have used the most naive definition, where we first 
tune the couplings of the potential V{\) to their values at the N = 00 multicriticality, and then take 
the N ^ 00 limit in the differential equation for the polynomials. What we implicitly have shown 
above is that the final result will depend on this precise procedure. This should come as no surprise, 
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since we in contrast to the conventional case of a finite p(0) really have to approach a critical point at 
some particular critical coupling(s). The whole analysis is then quite analogous to finite-size scaling 
in statistical mechanics [£^]- Finite N here plays the role of finite volume V, and as in statistical 
mechanics the notion of a (multi) critical point at finite A^ is ambiguous. One can define it to be 
precisely at the position of the infinite- A^ (multi)critical point, or one can choose other definitions (at 
the peak of the specific heat, or otherwise) that all agree in the N = oo limit. 

In order to remedy this situation, we clearly have to consider more general ways of reaching multicrit- 
icality. Instead of starting out with the couplings fixed precisely at the values corresponding to A^ = cxd 
multicriticality, we can consider a continuum of different schemes in which the couplings are tuned 
towards the multicritical values at just such a rate as to contribute non-trivial terms to the differential 
equation in the microscopic A^ — > oo limit. Because m couplings must be tuned to reach ?7ith-order 
multicriticality, this program is however not very practical from an analytical point of view. A much 
more economical procedure is to tune only one parameter, which in a generic manner takes us away 
from and toward the multicritical points. Following the conventions of matrix models applied to 2-d 
quantum gravity, we shall call this a double-scaling limit. It is, in a very precise sense, equivalent to 



finite-size scaling theory in statistical mechanics |2C] 



Suppose, then, that we have chosen our potential V{\) to be multicritical according to eq. ( p.l9| ) for 
some given M. Consider now the partition function defined by 

Z = fdMdet^^M e-W3)TrViM) _ (2.71) 



When g = 1 this coincides with our previous definition, eq. ( |2.1| ). We can tune g towards gc = ^ in 
such a way that 

y ^ Ar2--(g-i) (2.72) 

is kept fixed. This defines our double-scaling limit. The previous analysis was simply restricted to 
y = 0. Now all solutions will be parameterized by (functions of) this number y, which is kept fixed in 
the large- A^ limit. The appropriate scaling variable is then 

z ^ A^2™. 1^1 _ 2n^^j ^^.73) 

instead of eq. ( p. 45] ). This is completely analogous to the treatment of multicritical limits at the 



edge of the spectrum by Bowick and Brezin [27|. We can now probe a double-scaling region where 
g differs from unity only by an amount of the order _/v~^'"'^, while the eigenvalues A differ from zero 
by an amount of the order N~'^. Essentially all of the previous discussion goes through unchanged in 
this more general setting, with just some minor redefinitions. It follows from the above that the sole 
effect of including a non- vanishing y is to shift the scaling variable hy z ^ z + y. Accordingly, all of 
the previous discussion goes through unchanged in this more general setting, except that we should 
substitute the numerical values of the Painlev'e II transcendent at z = y instead of z = into the 
linear differential equation ( 2.38| ). 



2.7 A MULTICRITICAL MESOSCOPIC RANGE 

It is of some interest to compare the correct solution for the wave functions ipN^x) with what we would 
have obtained if we had ignored the apparently subleading terms in the defining differential equation 
(p. 38). Consider again for simplicity the case m = 1 and a = 0, for which we can compare the two 



solutions visually (see fig. 2). It is clear that although the two solutions differ in detail, especially very 
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close to the origin, their quahtative features are nevertheless quite similar. This has a simple algebraic 
explanation. The coefficient of the leading power of Ai\f{X) is proportional to k{m), which can be seen 
from eq. ( |2.19| ) and the definition (|2.9| ). Already for to = 1,2 k{m) is large compared to one, a fact 
which is improving for growing m: 

1 <C k{m) — > 4\/27rm2 for m ^ cc . (2.74) 

If now the leading power of A^^X) is dominating all the lower terms, we will get back the relation ( |2.16[ ) 
to the eigenvalue density as in the noncritical case and will therefore recover the "naive" differential 



equation from ( ^.28|) . The appropriate condition for x in this mesoscopic limit is 

1 < k{m) x^'" , (2.75) 

with the simplified differential equation reading 

Wn\x) il)j^{x)+ \cj^k{m)x H ^ \Vn{.x) = 0. (2.76) 

X \ X'^ j 

This mesoscopic differential equation can be solved analytically in terms of Bessel functions. Note 
that since it is an approximate equation it is no longer a-symmetric. At growing multicriticality m it 
will approximate the true differential equation ( |B.g| ) for smaller and smaller values of x, due to the 
behavior of k{m) eq. ( ^.74| ). 



As it has been mentioned eq. ( 2.76| ) can be mapped to the following equation of Bessel type 



where 

_ k{m) ,^_,, _ vrp(^-)(0) o^^, _ 2(-l)^a 

" 2(2m + l) "(2m + l)! ' ^" 2m + 1 " ^ > 

Boundary conditions are however not easily established. The solution satisfying the usual boundary 
conditions, namely regularity at the origin and normalizability, reads 



^n[x) ^ VXJ ^ (-i)iv(X) , (2.79) 



2m+l 



which still has to be normalized. (Since the differential equation cannot be trusted for small values of 
X, it is not obvious that the above boundary conditions are suitable; we shall for illustration stick to 
them here.) The approximate kernel and the approximate spectral density are then given by 

Kmeso{x,x') = Hui N'^d+i R^ ( N~2^ x, N~2^ x') (2.80) 

^ 2m+l ' 2 2m + l 2 2m + l 2 2m+l ' 2 J 



2(2m + l)! {x-x') 



and 
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(x) 



2(2m)! ] 2m + 1 



J 



_iW + J_ 



.W-J 



.(^)^. 



jW--/. 



2m + 1^2 2m + l 2 2m + l ^ 2 2m + l 2 2m + l 2 

Here the normalization has been fixed by requiring the matching condition 



^X)J.^ 



hm p^^™^ 



p 



(2m) 



(0) 



(2m)! 



k{m) 
27r 



(2.81) 



(2.82) 



The above "mesoscopic" spectral density of course has only approximate validity. For a = it simply 
coincides with the macroscopic spectral density near the origin (a feature it shares with the usual 
non-critical microscopic spectral density (|1.2D ), and is in this sense "exact". For a/{2m + 1) non-zero, 
but small, it does define an intermediate range where the set in of the universal mesoscopic oscillations 
is correctly encoded, and the macroscopic spectral density still not reached. To give an example of 
the approximation achieved, we show below, for a = l and m = l, the full microscopic spectral density 
(fig. 7) and the approximate mesoscopic spectral density (fig. 8). The full microscopic spectral density 
has been computed using the method described above, where, from the solution to the Painleve II 
equation, we find u+ = 2.141, u_ =4.218 and w = -10.087. 
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Figure 7: Microscopic spectral density for m=l criticality with a = \. 

r [x]/x2 

3r 




Figure 8: The approximate mesoscopic spectral density Pmeso{x) / x^ for m—\, a — 1. 
For m = (A;(0) = 4) the mesoscopic density reduces to the exact, noncritical microscopic density 
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3 Multicritical Chiral Unitary Ensembles 

Once the above results have been estabUshed, it is only a small step to extend them from the unitary 
ensemble to the chiral unitary ensemble. This is due to the fact that the partition functions, and 
their associated orthogonal polynomials, for the two different ensembles can be related to each other 
through appropriate shifts of the arguments. Let us recall |j|] that the pertinent matrix integral for 
theories in even dimensions, gauge groups SU{N>3), and fermions in the fundamental representations 
of these gauge groups, is described by complex matrices [^: 



dMdet^fM exp 
with an even potential 



— TrV{M^] 



(3.1) 



k 
Here M is a (2A''+|z^|) x (2A^+|i/|) block hermitian matrix of the form 



V{M^) = J2 T^^^ ■ (^-2) 




M = :,: (3.3) 



where W itself is a rectangular complex matrix of size A^ x (A^+|zv|). The translation of parameters into 
gauge theory language is that u represents the topological charge, and the space-time volume equals 
2A^+|i^|, the size of the matrix M. The integration measure in eq. ( p.lD is the Haar measure of W. 
This measure and the integrand of eq. (|3.1D both are invariant under left and right multiplications: 



W^TWU, U£U{N), T£U{N+\i^\). (3.4) 



It is convenient to rewrite the matrix integral ( |3.1| ) in terms of the eigenvalues of the hermitian matrix 
H = W^W 0: 

dWdW^det^f{W^W) exp\-NTrV{W^W)] , 

^ N 2 

n(dAiA?e-^^(^»))|deti,A}-i| , (3.5) 

(3.6) 



i=l 






detijzf-^ 



where a = Nf+\u\. Here the fact has been used, that one can switch from real positive eigenvalues Aj 
to real eigenvalues Zj P]. In the latter form, the matrix integral can be analyzed directly in terms of 
the unitary ensemble considered in the previous section. 

Let us introduce a set of even polynomials P^ {z'^) orthonormal with respect to the measure duj{z) = 
dz\z\'^°'~^^ exp(— A^y(z^)) on the real line. Then we can replace the Vandermonde determinant as usual, 
but this time by a determinant of only even polynomials -P^ (z^). In doing so we can directly use the 



even subsetP] of the orthonormal polynomials eq. (2^), which we have already calculated, by shifting 



a -^ a + ■^: 



Pi-\z') ^ Pir^)(.) . (3.7) 



We will also encounter the odd polynomials when calculating the kernel via the ChristofFel-Darboux formula. 
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In going from eigenvalues A to z in eq. ( |3.6D we have an even potential V{z'^) as we should, in contrast 
to the real positive picture (|3.5D. In order to justify the identification (|3.7D also in the multicritical 



case we still need a relation between the respective critical potentials. It is given by ||28|] 

Vra{X^) = 2y„(A) , (3.8) 

The critical potential for the chiral unitary ensemble has been calculated by imposing 

1 



Pmiz) = —k{m)z^^VT^^. (3.9) 

47r 



for the eigenvalue density on R, exactly as in the unitary case eq. ( 2. IS ). The factor of two in eq. 
( p.8| ) can be traced back to the difference in the saddle point equations of the two ensembles p^ . 
Consequently the wavefunctions of the chiral unitary ensemble at ?7i-th multicriticality can be taken 



from the unitary ones eq. (2.4) simply by shifting a ^ a + -' 



2- 



v;;°)(z2) ^ |zr+^e-fv'"^(^')p/")(z2) 

= |zr+^e-^^-(^)4"+^^(z) 



4r^^(z) , (3.10) 



where we have identified the couplings gk = gu according to eq. (|3.q ). The wavefunctions of the 

unitary ensemble "02; ^ i^) '^o'^ have to be evaluated at 2N instead of A''. They are determined by 
the corresponding universal differential equation as being discussed in great detail in the previous 
section. 

The eigenvalue density in the microscopic limit is obtained from the kernel expressed in terms of even 
and odd wave functions ipk ^ i^)'- 

N-l 

Km{z\w^) = |zu;|°+ie-i^(^'-(-')+^'-("'')) ^p/")(z2)p/°)(u;2) 

z=o 

..(°+^)/.A,,.("+5)/ 



Y.^r~'\^w2' ^i^) 



^21 
1=0 



= C2NC2N~1 ^2 _ ^2 [^2N {z)V2N-2 M " V2N-2 iz)i^2N (^) 

Here we have made use of the Christoffel-Darboux identity and the recursion relation for the wave 

functions. The advantage of the last expression is that we have avoided to expand '4'2N-2 i-^) ^^ 1/A^. 
Imposing the same scaling for the real eigenvalues z as in eq. ( |2.31 ), 



1 



N^^^^z (3.12) 



we obtain 



Ps(x) = lim lim N ^m+iK^ (x'^N ■^m+i^y^N ^m+A (3.13) 

N^oo Z \ X 
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3.1 The multicritical mesoscopic range 



Since in the unitary case we could obtain analytic expressions for the wave functions in an inter- 
mediate, mesoscopic range we can do the same here, using the results of subsection ^J. The range 
of applicability is thus given again by eq. ( 2.75| ). The unnormalized, analytic expressions for the 
mesoscopic wave functions eq. ( ^.791) , which are needed here read (see our discussion about boundary 
conditions in the previous section): 



where 



X = 



k{m) 



^2m+l 



27V ^ 
2N-1 ^ 



27rp(2™)(o) 






^2m+l 



a 



2a + 1 
2(2m + l) 



(3.14) 



(3.15) 



2(2m + l) ""'' '"'(2m + l)! 

Inserting these expressions into eqs. ( |3.11 ) and ( 3.14| ) we obtain the following expressions for the 
approximate mesoscopic kernel and spectral density 



Kmeso{x,x') = lim N 2"^+'^Kn (x'^N 2m+i^a;'^Af 2,«+i) 

Af— +00 V / 



m+4 



27rp(^'")(0)(xrr^: 
(2m + 1)! a;2-x'2 



(x J^_i(X')Vi(X) - x' J^_i(X)J^^i(XO) (3.16) 



and 



_ „(2m,) (r\\rf,2m 
Pmesoix) = \,:^Z ( ^ ( Vi (^)' + "^/J-i (^)' 



2 (2m)! 
The normalization has been fixed by requiring again 



hm p^^™^ 



-i(2m) 



(0) 



2n 

^J.,i(X)J. i(X) 

k{m) 



(2m)! 



47r 



(3.17) 



(3.18) 



For m = the above mesoscopic density ( 3.17| ) reduces to the exact noncritical microscopic density 

0. 



In fig. 9 we exhibit the exact microscopic density (3.14) for m = 1 and a = 0, constructed by solving 
the Painleve II equation for a = 1/2, and the mesoscopic approximation to it eq. ( 3.17] ). 
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Figure 9: Microscopic spectral density for the ra — \ critical chiral unitary ensemble with a — Q (bold line) 
and its mesoscopic approximation (real line). The normalization is determined by requiring Imix^oo p'si^) — 
p"(0) =8/7r (dotted line). 
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4 Conclusions 

In this paper we have investigated random matrix models which may be related to QCD-like gauge 
theories with unbroken chiral or flavor symmetries. From the Banks-Casher relation |y] the vanishing 
of the corresponding order parameter ('(/'V) implies a vanishing spectral density of the Dirac operator 
at the origin. The case where this density vanishes as an even power can be mapped to the multicritical 
points in the appropriate matrix models if one follows the same arguments as in |jll for the broken 
phase. 

The new universal results we find for multicritical matrix models in the microscopic large- A^ limit may 
suggest a universal behavior for the microscopic density of the Dirac operator itself. In contrast to the 
case of broken symmetries, there are here no known spectral sum rules of these theories with which to 
compare the results from random matrix theory. The relevance of these results for the Dirac spectrum 
of theories with unbroken chiral/flavor symmetries is therefore still far from obvious, and we can at 
present only view it as an intriguing conjecture. It would be particularly interesting to investigate if 
an approach analogous to that of a chiral Lagrangian [^ exists in this case. 

The detailed universality predictions, the listing of universality classes, and the "critical exponents" 
V = l/(2m + 1) can all be tested in lattice gauge theory with an appropriate fermion prescription. In 
contrast to the noncritical case |^, the universal results for the microscopic spectral densities are no 
longer given in terms of elementary functions. Instead, we derive universal differential equations for 
the orthogonal polynomials of the multicritical models, which then yield the spectral density from its 
kernel. For the practical purpose of comparison with lattice gauge theory data, the generalization of 
the present formulation to that of massive fermions may be highly advantageous. Recently the massive 
case for non-critical potentials was worked out in detail for the chiral unitary ensemble and proven to 



be universal [^ in the appropriate scaling limit (see also ref. |3C] for an explicit computation in the 
Gaussian case) . The same universality of the massive case extends also to unitary ensemble [^] . The 
multicritical cases for both ensembles with non-zero masses can be worked out completely analogously 
using the wave functions derived in this paper. The universality arguments also generalize accordingly. 

We have discussed here in great detail the solution for first order multicriticality with an arbitrary 
number of massless fermions, giving evidence for the universality of the determining differential equa- 
tion and thus for the microscopic spectral density. For higher order criticality we have derived the 
general form of the differential equation and conjectured its universality. With growing multicritical- 
ity m the analysis becomes increasingly involved, but all the necessary ingredients for carrying the 
computations through for arbitrary m have been presented here. Moreover, in a certain mesoscopic 
large- A^ limit we have been able to give approximate, analytic expressions for the universality classes 
of the orthogonal polynomials as well as the microscopic density for any order of multicriticality and 
any number of fermions. 

In the noncritical case, the relevant scale that determines the microscopic limit is the spectral density 
at the origin, /o(0). All eigenvalues are in that case blown up near A = by the scale p{0)N. As could 
have been guessed from the outset, in the mth multicritical case the relevant scale is set by the first 
non- vanishing derivative of the macroscopic spectral density at the origin, p^'^"^^{0). All eigenvalues 
must be blown up by a scale determined by this number. In addition, and this is what requires a 
more detailed analysis, we have seen that the eigenvalues must be blown up with the volume N raised 
to the exponent u. The noncritical case is simply that of ?7i = 0. Viewed from this perspective, the 
universality classes derived in this paper are the completely natural generalizations of those derived 
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previously ^ . Universality simply translates into the freedom of redefining only one scale, p(2'^)(0), 
in order to leave all results invariant. 



One might think that the most obvious application of the universal results described in this paper 
would be in connection with chiral symmetry restoration at finite temperature or density. Because one 
parameter needs to be tuned in each case (temperature T, or chemical potential fi), one could believe 
that these two cases, when occurring at continuous phase transitions, would be good candidates for 
the m = l universality class. We are mildly skeptical about this, and have instead, in the introduction, 
suggested what seem to us to be more promising candidates. There has already been a large amount 
of work on possible extensions of the by now essentially established T = results to the T ^ and 
fi ^ situations (see, e.g., ref. for just a limited selection of papers, and also others cited in ref. 
11). 
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A A symmetry among the orthonormal polynomials 

In this appendix we prove a lemma relating polynomials orthonormal with respect to measures dXw{X) 
and dXX'^w{X). It follows from a more general statement in the book of Szego ([^, theorem 3.1.4). 
We have referred to this relation as the a-symmetry in the main text. 

Lemma: Let Pjl"^(A) denote the polynomials orthonormal w.r.t. the measure 

(iAu;(")(A) = dA|A|2°exp(-W(A)) , a > , AgR, (A.l) 

where V{X'^) is an even polynomial with positive highest coefficient. The following identity then holds: 

(A) 



A 



pr^'w ■ (A.2) 



~( +1) p(") 

Proof: Since (A) = ^"^^ — is an even polynomial of degree 2n it is orthonormal to all odd 

polynomials w.r.t. dXw^"~^^\X) . Because of 

dXw'^"+^\X)Pt^'\x)Pt^'\x) = / dX ^;(^+i)(A) ^2'»+i^^ pn+iiAj 



oo 



A2 

dX t.(")(A)PS+i(A)p£|,(A) 



one has 

^^"'"'^(A) = Pt^'\x) . (A.4) 
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Inserting the above result into the definition (|2.4| ) of the wave functions '0n(A) immediately leads to 
the a-symmetry as stated in eq. (2.30). 



B The m > 2 multicritical case 



In this appendix we give the explicit result for the asymptotic form of the differential equation ( 2.13 ) 
for m = 2, and conjecture the form for general m. The m = 2 minimal multicritical potential reads 



V2(A) 



1 



1 



1 



6^ 



with 



92 



1 

2^ 



94. 



16 



96 



32 , 



(B.l) 



(B.2) 



Performing the calculations exactly along the same lines as in subsection 2.3 we get the following exact 
expressions 

An{X) = NCn {92 + 54(4+1 + cl) + 56 + + C^) + C^(c^+1 + + cLl))) 

+Ncn{gi + 9(i{cl+i + cl))X^ + NCn9&X^ 
Bn{X) = Ncl[g, + g^{cl_,^ + cl + cl_{))x + Nclg.X' + {l-{-ir)j (B.3) 

Inserting them into the exact from for Gn{X) eq. ( p. 15 ) and rescaling variables according to 

X = mX , (B.4) 
we expect that the following non-trivial contributions to be of order 1: 

V = Ni {g2 + 2g4C% + 2gQC%{c%^^ + + 

p = m[g4 + 2gec%) (B.5) 



This can be checked explicitly from the string equation for V2(A) when inserting the ansatz ( 2.44 ) for 
m = 2. We do not display the differential equations for the scaling functions /2(-z), 52 (-2), h2{z) etc. 
here. Finally the asymptotic expression for the differential equation of the wave functions reads in the 
microscopic large- limit 



with 

N-^/'>Gn{x) 



^'^{x) - N-^"'Fn{x)^'n{x) + N-^''>Gn{x)iI>n{x) 



2q^x + 'igQX'^ 







(B.6) 



u+ + q+x'^ + gex^ 



3 1 I \ 

c%u+u^ + ((-l)^a - -)v + icUu+q^ + u^q+) - — + {{-if a - -)p\ x^ 

{2u+ + q+x^){v + px^) + (-l)^2a(g+ + 2gQX^) , 3 2 8 , {-ifa-a^ 
+ : — A \-CNgQX H ^ (B.7j 



u+ + g+x^ + gex^ 



x^ 
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where only the dominant contributions are displayed. Moreover we have used that the quartic and 
sextic terms in Gm{x) are subdominant due to the following identities, which are valid up to subleading 
terms: 

= 2c%gQ{u+ + U-.) + 2c%q+q- - pv 

= 4c^(76(g+ + 9-)-/ . (B.8) 



The former follows directly from the string equation (|2.7D whereas the latter can be seen to be valid 
when expressing it in terms of v. 

For the general case of mth multicriticality we obtain the following form: 

iP';f{x)-N-^FN{x)ilj'j^{x)+N~^GN{x)i'N{x) = 0. (B.9) 
The -F/v-term does not simplify: 

2 

1 /V 2m + l A\r( t) 

N-WFn{x) = , ^^"^^ , (B.IO) 

AT 2m + lAN{x) 

1 1 

where ^"''+^ Aj^{x) is a polynomial in x = A^2m+i^ of the order 2m. The GTv-term which we 
conjecture to be of the following form will simplify considerably: 

N-^^Gn{x) = Q{X) + —-^ + C^52„^+2X^'" + • (B.ll) 

jSf 2m+iAN{x) ^ 

Here Q{x) and R{x) are both polynomials of order 2m — 2 in x. They depend on certain universal 
combinations of coupling constants and recursion coefficients, as it has been explicitly shown in the 
two examples for m = 1 and m = 2 in eq. (I^D and (|B7| ). The latter can be brought into the 
conjectured form when dividing out the quartic term in the denominator. 
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